Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 5 February 2008 (MN WF$t style file v2.2) 



Disk Galaxy Evolution Along the Hubble Sequence 



Marios Kampakoglou and Joseph Silk 

Department of Physics, University of Oxford, Keble Road, Oxford OX1 3RH, 



United Kingdom (e-mail: mariosk, silk@astro.ox.ac.uk) 



o 
o 



CO 



CN 
> 

o 

\o 
o 

o 

0,; 

6 

CO 
C3 



13 



5 February 2008 



ABSTRACT 

Galaxy disks are characterised by star formation histories that vary systematically 
along the Hubble sequence. We study global star formation, incorporating supernova 
feedback, gas accretion and enriched outflows in disks modelled by a multiphase in- 
terstellar medium in a fixed gravitational potential. The star formation histories, gas 
distributions and chemical evolution can be explained in a simple sequence of models 
which are primarily regulated by the cold gas accretion history. 



1 INTRODUCTION 



The two most striking characteristics that define the Hubble 
sequence are morphology and star formation activity. The 
former is best addressed by numerical simulations (Abadi et 
al. (2003); Sharma & Steinmetz (2005); Robertson et al. 
(2004); Governato et al. (2006)). However the latter aspects 
are so complex that most discussions of disk star formation 
and chemical evolution are based on analytical calculations 
(Efstathiou (2000); Silk (2001); Ferreras & Silk (2001); 
Silk (2003); Matteucci (2006); Naab & Ostriker (2006)). 
This paper extends the analytic approach to study star for- 
mation histories that vary systematically along the Hubble 
sequence. We study global star formation, incorporating su- 
pernova feedback, gas accretion and enriched outflows in 
disks modelled by a multi-phase interstellar medium in a 
fixed gravitational potential. 

One of the current problems afflicting galaxy formation 
models is the role of gas infall. For both ellipticals (Bower 
et al. (2006); Croton et al. (2006)) and massive disks at 
z ~ 2 (Forster Schreiber et al. (2006)), infall rates are suf- 
ficiently high from CDM theory (ellipticals and disks) and 
observations (disks) that infall must be quenched, otherwise 
distant ellipticals are too blue and the disks are too mas- 
sive by the current epoch. The problems may be related via 
feedback from AGN, but the details are poorly understood 
with regard to the resulting star formation history, gas infall 
rates and chemical evolution. In this paper, we focus on disk 
galaxies, and develop a phenomenological description of disk 
evolution in which the onset and duration of the gas infall 
history are found to be the controlling parameters. The star 
formation histories, gas distributions and chemical evolution 
can be explained in a simple sequence of models which are 
primarily regulated by the cold gas accretion history. 

One of the main features of galaxies that characterises 
the significance of the Hubble sequence is the wide range 
in young stellar content and star formation activity. This 
variation in stellar content is part of the basis of the Hubble 
classification itself, and understanding its physical nature 
and origins is fundamental to understanding galaxy evolu- 
tion in its broader context. In this paper, we construct a 



sequence of evolutionary models in which present-day prop- 
erties of different types of disk galaxies are reproduced. 

The first implementation of supernova-driven feedback 
in the context of CDM was to account for the properties of 
dwarf galaxies (Dekel & Silk (1986)), and subsequent stud- 
ies have explored the role of supernova feedback in more 
massive galaxies. Feedback is an important element in our 
attempts to model galaxy evolution. Energy injection from 
supernovae is probably the most plausible feedback mecha- 
nism for systems with virial temperatures higher than 10 5 
K. Winds from quasars might also disrupt galaxy formation 
or limit the growth of central black holes. Here, we will be 
concerned exclusively with supernova-driven feedback and 
we will not consider feedback from an active nucleus. The 
model presented in this paper is similar to the simple self- 
regulating model with inflow and outflow developed in Efs- 
tathiou (2000), the main difference being in the infall model 
that we implement. We adopt an exponentially decreasing 
infall rate normalised in order to reproduce the observed to- 
tal disk mass density in the solar neighbourhood, whereas 
in Efstathiou's model it is the conservation of specific an- 
gular momentum that specifies the final radius in the disc 
for each gas element. Angular momentum conservation is 
known to be a poor approximation, at least in the numer- 
ical simulations, and our more phenomenological model is 
easily adapted to the chemical evolution constraints. 

The main result of this paper is that star formation 
histories of the different types of disk galaxies can be re- 
produced using a one parameter model. The key parameter 
of our model is the time corresponding to the onset of the 
infall. Using this free parameter, we reproduce the distribu- 
tion of disk birthrate parameters b, the ratio of the current 
SFR to the average past SFR, for each type of disk galaxy 
as presented in Kennicutt et al. (1994) and summarised in 
Figure [1] We then explore the various implications of the 
model for the radial and temporal dependences of the gas 
fraction, star formation rate and metallicity. 

The layout of this paper is as follows. Section briefly 
reviews the main points of our model and presents some ba- 
sic results. The model is extended in Section [3] to include 
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an improved treatment of winds and chemical evolution. Fi- 
nally, in Section [4] we discuss our results and present our 
conclusions. 



2 MODEL 

The model described here is a self-regulating model with 
inflow and outflow. The disk is considered to be a system 
of independent rings each of 35 xr^ pc wide (where r& is 
in units of kpcj]]. Neither radial inflows nor radial outflows 
are considered. The ring centred at the galactocentric dis- 
tance 7-q = 8.5 kpc is labelled as the solar neighbourhood. 
For the present-day total disk surface density in the solar 
neighbourhood, we adopt a value E tot (rQ, t g ) = 60 M pc -2 
(Holmberg & Flynn (2004) found 56±6 M Q pc" 2 for their 
disk model) . 



2.1 Dark Halo and Disk Model 

The dark halo is assumed to be described by the Navarro 
Frenk & White (1995) profile, 

Sp c 



p(r) 



(Cx)(l + Cx) 2 



(1) 



where x=r/r v , p c is the critical density, r v is the virial radius 
at which the halo has mean overdensity of 200 with respect 
to the background and C is the concentration parameter. 
For the present model, we adopt a value for C of 10. The 
circular speed corresponding to this profile is 



vh (r) 



2 1 [ln(l + Cx) - Cx/(l + Cx)] 
" v x [ln(l + C) - (7/(1 + C)] 



(2) 



where v^,= G ^ Iv , and M v is the mass of the halo within the 
virial radius. 

For the surface density of the disk at t = 0, we assume 
an exponential profile: 



E tot (r,0) = E e~ 



(3) 



where Md = 27rr- 2 Eo, Et ot (r, 0) is the total surface density 
of the ring centred at galactocentric radius r at t= 0, Md 
is the total disk mass at t= and is the scale length. So 
in our model the entire gas disk has formed instantaneously 
at t = 0. The rotation curve of a cold exponential disk is 
given by (Freeman (1970)): 



v D (r) = 2v c y [I (y)K (y) - h(y)Ki(y)] 



(4) 



where w=i— and vi= GM ° . For ratio — we adopt a value 
equal to 0.45 that corresponds to the Milky Way. Models us- 
ing standard disk formation theory with adiabatic contrac- 
tion within the cuspy halo (Klypin et al. 2002) reproduce 
the broad range of observational data available for the Milky 
Way by adopting a virial mass equal to M„ » 1O 12 M0, bary- 
onic mass equal to Mb ar « 4 — 6x 1O 1O M0 and virial radius 
r v — 258 kpc. These values in combination with an adopted 
value for the disk scale length equal to r^ = 3 kpc (Sack- 
ett (1997)) give |* = 0.45. In this work, we keep the ratio 



1 For a disk galaxy with scale length 
each zone is equal to 35 X 3 = 105 pc 



3kpc the width of 



— constant for all the models we present PI. The ratio of 

v c 

virial radius of the halo to disk scale length is defined as the 
collapse factor f C oll=^ = 50; such a value is needed to re- 
produce a median value of « 0.05 for the dimensionless spin 
parameter of the halo \h (Efstathiou (2000)). The disk is 
truncated at radius r/r^ m 7. 



2.1.1 Two-component(stellar and gas) rotating disk 

The stellar radial velocity dispersion a* is related to that of 
the gas clouds: 



(J* = CWT 



(•>) 



Here we assume a — 5 (Efstathiou (2000)). The scale- 
height for a two-component rotating disk is (Talbot & Ar- 
nett (1975)): 



H = 



ttGE, (1 + EiTEi) 



(6) 



2.2 Epicyclic frequency Model 

The epicyclic frequency is given by the expression 

„ 1 r did . i In ,_ v 

k= Ml+2-^) 1/2 , (7) 

in which ui is the angular velocity of the disk. In the above 
equation we replace lo — v to t/r where n 2 ot = v% + v D . So 
the epicyclic frequency(in units of 10 _15 sec _1 ) is: 

K = 0.035^2^(1 + (8) 

r vtot dr 



2.3 Star Formation Rate(SFR) 

In this work we adopt the star formation law proposed in 
Wang & Silk (1994) , to which we refer for a detailed descrip- 
tion. The star formation rate (in units of Mopc -2 Gyr -1 ) is 
given by : 



1>(r,t) = 



6kE 9 (1-Q 

Q 



2U/2 



(9) 



where k is the epicyclic frequency, E 9 is the gas surface 
density in units of Mope -2 , Q is the gravitational instability 
parameter(Toomre (1964)) and e is the efficiency of star 
formation. For the model in this paper we set e = 0.02. 

The star formation rate defined above has the follow- 
ing features: (i) star formation can only occur if the disk is 
gravitationally unstable, Q < 1; (ii) the degree of instability 
of the disk, measured by Q, is directly linked to the star- 
formation rate; the smaller Q is, the more rapidly stars are 
formed. 

The stability criterion for a two-component (stellar and 
gas) rotating disk can be written as (Wang & Silk (1994)): 



7TG <J„ (7* 



(10) 



The ratio defines the disk scale length r^. The disk scale 

C I 3«D Nl/3 



length is r d 



4x200p c / coii (^.)2 
fcoll 



kpc 
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where E* and <r* are the stellar surface density (in units 
of M pc~ 2 ) and the radial velocity dispersion (in units of 
km sec -1 ), and k is the epicyclic frequency. 

2.4 IMF 

The adopted stellar initial mass function is of the standard 
Salpeter form : 

dN-k 



= Bm mi < m < m u , x — 1.35 



(11) 



dm 

where mi — 0.1 M ,m u = 50Mq. 

For the IMF adopted, one supernova is formed for every 
125Mq of star formation, assuming that each star of mass 
greater than 8M releases 10 -E51 ergs in kinetic energy in 
a supernova explosion. Therefore, the energy injection rate 
per unit surface area (in units of erg s~ 1 pc -2 ) is given by: 



2.5xl0 32 £ 5 ie c <!/> 



(12) 



where the parameter e c defines the percentage of supernovae 
energy that goes to the ambient medium. For this model we 
adopt e c = 0.03. 



2.5 Energy dissipation due to cloud inelastic 
collisions 

The rate of energy loss per unit surface area (in units of 
ergs^pc -2 ) is given by (Efstathiou (2000)): 



5.xl0 29 (l + 



a* /a 9 



)o- 5g E 



5g- 



(13) 



Using the above equations we know the energy balance in 
each ring of the disk. We define the radial velocity dispersion 
of clouds (in units of kms -1 ) as: 



a g {r,t) = 



2xE(r,t) 
M g (r,t) 



(14) 



where E(r,t) is the energy balance at galactocentric radius 
r and time t, and M g (r,t) is the mass of the cold gas at 
galactocentric radius r and time t. 

In summary, the energy balance in each ring of the 
disk defines the cloud and stellar radial velocity dispersion 
a g ,<j* through equation (14) and equation (5). Combining 
this piece of information with the cloud and stellar surface 
densities E 9 ,E*, we calculate the degree of instability at 
galactocentric radius r through equation(lO). The star for- 
mation rate is regulated by the degree of instability, so the 
final form of the star formation rate is derived by inserting 
equation (10) into equation (9). 

2.6 Infall Model 

The adopted form for the gas infall rate is an exponentially 
decreasing function in which the rate of gas infall (in units 
of Mopc -2 Gyr -1 ) in each ring is expressed as: 

f(r,t) = A{r)e^ (15) 

where tj (in units of Gyr) is the infall time scale. The in- 
fall rate f(r, t) is normalised to the present-day local disk 
density, f* 9 f(r,t)dt= E tot (r, t g ), where E tot (r, t g ) is the 
present-day total disk surface density of the ring centred at 



galactocentric radius r, and t g is the age of the galactic disk 
(t g = 10 Gyr). Assuming that the present-day total masses 
of the different rings exponentially decrease with increas- 
ing galactocentric distance with scale- length r<j, the form of 
A(r) can be written as (Chang et al. (2002)): 



A{r) = 



T, to t{rQ,t g ) - E tot (7-0,0) 



rt g 

J t,. 



(16) 



T J dt 



where E tot (r0, 0) is the total disk surface density in the solar 
neighbourhood at the epoch which the formation of the disk 
begins, and E tot (r0, t g ) is the present-day total disk surface 
density in the solar neighbourhood. We adopt an exponen- 
tially decreasing infall rate normalised in order to reproduce 
the observed total mass density in the solar neighbourhood. 
In our infall model, we have introduced a parameter ti ow 
that corresponds to the time that infall switches on. Ad- 
justing this parameter we can reproduce the distribution of 
disk birthrate parameter b, the ratio of the current SFR to 
the average past SFR, for each type of galaxy presented in 
Kennicutt et al. (1994). The infall time scale 77 (in units 
of Gyr) is also a free parameter of the model but we keep it 
constant, adopting a value of 0.5 Gyr. 

According to the adopted infall form, since 17 does not 
vary with radius, the disk keeps the exponential form with 
scale length r d at each epoch of evolutionary time. 



2.7 Hot Phase Model 

An expanding supernova remnant will evaporate a mass of 
cold gas. The rate of the evaporated mass per unit area is 
(Efstathiou 2000): 



M"^lxl0~ 7 - 



E 59 (l + ^) 



. . r>0.71 0.29 t-,0.71 /.-0.29 

■x5i3 7 #5i /e 



(17) 



in units of Mepc^yr -1 , where Es 9 is the gas surface density 
in units of 5 Mqpc -2 , (T5 9 is the cloud radial velocity disper- 
sion in units of 5 kms -1 and S13 is the supernovae rate in 
units of 10 _13 pc _3 yr _1 . The parameter 7 relates the blast 
wave velocity to the isothermal sound speed (v& = jch, 
7 w 2.5). The parameter /e is defined through the evapora- 
tion parameter (in units of pc 2 ): 



280 



0"5g 



■ J 5 g 



= /eS £ 



(18) 



(1 + if) 



where E ot q is the evaporation parameter in the local solar 
neighbourhood and <f> K is a parameter that quantifies the 
effectiveness of classical thermal conductivity. For this model 
we adopt E elI Q = 95 pc 2 and 4> K = 0.1. For a system with 
porosity Q close to unity, the temperature and the age of a 
supernovae remnant are given by: 



to = 5.5xl0 D Si3~ 



5/11-6/11 fj.-3/ll 3/11 



7 



£ 5 7 



yr, 



T = 1.2x10* 5*^7 



4 ^6/11^-6/11^8/11^-8/11 



L K, 



The density of the ambient hot phase is 



A Ov, l n-3 cO-36 -0.36 T7.0.61 /•- 0.393 „-3 

n h = 4.3x10 Si 3 7 E 51 / E cm . 



(19) 
(20) 

(21) 
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For 10 5 <TX10 6 .ft r we adopt a cooling rate of A 



2.5xlO" 22 T 



ergcm s 



So for a gas with primordial 



composition, the cooling time t coo i is 



= 2t„n 



1 fO.t 



(22) 



2.8 Outflow Model 

As mentioned in the previous section, supernova bubbles 
expand, evaporate cold gas and compress the ambient in- 
terstellar medium(ISM). The compressed ISM will also be 
driven by a form of wind. In order to study the effects of 
outflow, we assume a simple phenomenological model for 
galactic winds where the wind mass-loss rate Mw (in units 
of M0pc _2 Gyr _1 ) is assumed to be proportional to the star 
formation rate. 



Mv 



kip 



(23) 



We explore many values of k between and 1. The hot 
gas that escapes from the halo is removed permanently. Su- 
perwinds indeed are found around Milky Way-type galaxies 
(Strickland et al. 2007). This is puzzling from the theo- 
retical perspective because supernovae are considered to be 
incapable of driving a strong wind from the gravitational 
potential well of a massive galaxy. Also, at least one exam- 
ple is known of an extended x-ray halo around a normal 
massive spiral galaxy (Pedersen et al. 2006). In this case, 
gravitational accretion and shock heating provides the most 
plausible source for heating the gas. 



2.9 Galactic Fountain Model 

Models of gas flow on the galactic scale were first intro- 
duced by Shapiro & Field (1976) and subsequently devel- 
oped by Bregman (1980) and others. The galactic fountain 
originates from the supernovae that warms up the disk gas 
to temperatures of 10 6 K. The upflowing gas cools and con- 
denses into neutral hydrogen clouds that rains onto the disk. 
The models assume the height to which the hot gas will rise 
and the expected rate of condensation in the cooling gas de- 
pends only on the temperature of the gas at the base of the 
fountain and the rate of cooling of the upflowing gas. In the 
model described in this paper the hot gas that is not lost 
from the disk returns to the disk at the radius from which 
it was expelled after a time t coo i. 



Input parameters 



Solar radius 

Total disk surface density at r@ 

Star Formation IMF 

Solar metallicity 

Effective yield 

Star formation efficiency 

Infall timescale 

Galactic disk age 

Thermal Conductivity effectiveness 
Solar evaporation parameter 
Disk to Halo velocity ratio 
Collapse factor 
Wind Model Parameter 



re = 
^tat(r@,t g ) 



8.5 kpc 
= 60 M Q pc" 



Salpctcr 
Z Q = 0.02 
p = 0.02 
e = 0.02 
T f = 0.5 Gyr 
10 Gyr 
= 0.1 
= 95 pc 2 
= 0.45 
= 50 
0.2 



to 



fcoll 

k = 



Table 1. Input parameters in common for the five models we 
examine in this paper. 



Type 


Disk Mass M d (Mq) 


r d (kpc) 


Sa 


5.5 x 10 1(1 


3. 


Sb 


4. x 10 10 


2.7 


Sc 


2.7 x 10 10 


2.37 


Sd 


1.5 x 10 10 


1.95 


Sm 


1. x 10 10 


1.7 



Table 2. Five models adopted in this paper 



2.11 Different Types of Spiral Galaxies 

To explore the evolution of different types of spiral galaxies 
we adopt five models in this paper (See Table [2]). These 
models differ only in the mass of the disk and the disk scale 
length Td- All the input parameters that these models have 
in common are summarized in Table [1] 



2.12 Basic Results 

In this section will present the basic results of the model dis- 
cussed in the previous sections. Unless otherwise stated, we 
adopt a value 0.2 for the parameter k in the outflow model. 
For clarification of the plots, when results from models Sc 
and Sd are similar we define a new type of disk galaxy Sc/d. 
The data for Sc/d type come from the average between the 
Sc model data set and the Sd model data set. 



2.10 Chemical Evolution 

We include chemical evolution in the model using the in- 
stantaneous recycling approximation. So we assume that all 
processes involving stellar evolution, nucleosynthesis and re- 
cycling take place instantaneously on the time scale of galac- 
tic evolution. The equation of galactic chemical evolution is 

E g dZ = pipdt + (Z F - Z)f (24) 

where f is the infall rate (in units of MQpc _2 Gyr _1 ) and p 
is the yield (Pagel (1997)). We adopt a yield of p = 0.02 
and assume that the mass accreted to the disk has zero 
metallicity (Zf = 0).We normalise the metallicities to the 
solar value for which we adopt Zq = 0.02. 



2.12.1 Distribution of disk birth rate b 

Values of the birth parameter b were calculated for each type 
of spiral galaxy. For the calculation of b values we ignore the 
SFR that corresponds to the first 10 8 yr, due to a deficiency 
in our model. More specifically in the model we assume that 
the entire gas disk has formed instantaneously at t = 0: 
this is unrealistic and leads to very high rates of star for- 
mation for the first 10 8 yr. Figure [5] shows the disk birthrate 
parameter b subdivided by galaxy type as a function of our 
free parameter ti ow that corresponds to the time that infall 
switches on. The dotted line indicates the median value of 
birthrate parameter while the dashed lines show the quar- 
tiles of birthrate parameter presented in Kennicutt et al. 
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Type 


b(median) 


tiow med (Gyr) 


Sa 


0.07 


0.10 


Sb 


0.37 


7.00 


Sc 


1.04 


8.05 


Sd 


0.63 


7.80 


Sin 


1.31 


8.35 



Table 3. Median value of b and t[ ow value for different models. 



E 




■low" Mow_med 



Figure 1. Gaussian fits to the b parameter dispersion for each 
type of disk galaxy 



(1994). In Table O we present the value of the free parame- 
ter ti om we adopt in order to reproduce the median value of 
birthrate parameter for each type of disk galaxy. 

To make our results more useful for observers, we fit 
Gaussians to the 6 parameter dispersion for each type of 
disk galaxy (Figure (TJ. The full width at half maximum for 
each different Gaussian equals the width of the ti ow param- 
eter distribution in order to reproduce the quartiles of the b 
parameter distribution presented in Kennicutt et al. (1994). 
In Table |4j we give the value of a for the Gaussian fits as 
well as the width in ti ow in order to reproduce the quartiles 
of the b parameter distribution. The width of the b param- 
eter distrbution between the quartiles for each type of disk 
galaxy is presented in the last column. 



Type a At low (Gyr) Ab 



Sa 


1.716 


4.04 


0.1 


Sb 


0.386 


0.91 


0.3 


Sc 


0.205 


0.48 


0.6 


Sd 


0.183 


0.43 


0.4 


Sm 


0.270 


0.64 


2.0 



Table 4. The value of a for the Gaussian fits (second column), 
the width in t[ ow in order to reproduce the quartiles of the 6 
parameter distribution (third column), the width of the b param- 
eter distribution between the quartiles of the distributions (fourth 
column). The values are shown for different types of disk galaxies. 



2.12.2 Density Profiles 

Figure [3] and figure [4] show the evolution of the gas and 
stellar surface densities, respectively, (in units of MqPC -2 ) 
for different ages and different type of disk galaxies. The 
results are shown for ages of 0, 0.1, 1, 3, 6 (top - bottom) and 
10 Gyr [dashed line). The star formation rates are initially 
high and hence the time scale for star formation is short. 
Almost the total amount of gas in the central region of disk 
is transformed to stars in 10 8 yr. The star formation rate 
declines rapidly after 1 Gyr. As figure [3] and figure [4] show 
the star formation at early times is concentrated to the inner 
parts of the disk which have a high surface density, and hence 
the gas distribution develops a surface density profile with 
an inner 'hole', similar to what is seen in the HI distributions 
in real galaxies (Deul & den Hartog (1990)). 

Model Sa has parameters similar to those of the Milky 
Way. Gas surface densities from direct observations at the 
solar neighbourhood is 8±5 Mq/pc 2 (Dame (1993)) in 
good agreement with the predictions of the model. Slightly 
higher values of ~ 13 — 14 M©/pc 2 have been reported by 
Oiling & Merrifield (2001) and are in good agreement with 
Sb model. Taking into account the scatter in the measure- 
ments of Milky Way scalelength (rd ~ 2.5 — 3.5kpc Sackett 
(1997) ) Sb model can also be a good candidate for the Milky 
Way. 

Direct observations of the stars at the solar neighbour- 
hood of Milky Way by Gould et al. (1996) and Zheng et 
al. (2001) using HST observations of M stars in the Galac- 
tic disk determine a surface density of 12.2 — 14.3 M©/pc 2 . 
Visible stars other than M stars contribute ~ 15 Mq/pc 2 , 
resulting in a total stellar surface density in the range of 
~ 27 — 30 Mq/pc 2 . This value is in good agreement with 
predictions from Sa and Sb models. 



2.12.3 Gas Velocity Profiles 

Figure [S] shows the radial distribution of the gas velocity 
at 10 Gyr for different types of disk galaxies, Sa (solid 
line), Sb (thick solid line), Sc (dashed line), Sd (thick dashed 
line) and Sm (dotted line). Gas velocity profiles are in 
good agreement with recent theoretical predictions (Ri- 
cotti & Ferrara 2002) that in a multiphase low metallic- 
ity (Z«5 x 10 _3 Zq) interstellar medium the velocity prob- 
ability distribution should be close to a Maxwellian with 
velocity dispersion cr^ 1 lkm/sec - 1 . Furthermore studies of 
nearby face-on galaxies show that the velocity dispersion in 
the HI layer is decreasing monotonically from about 10 — 13 
kms -1 in the optically bright inner regions to 6 — 8 kms -1 
in the very outer parts (Kamphuis & Sancisi (1993)). 



2.12.4 Star formation Profiles 

Figure [6] shows the evolution of the radial distribution of 
the star formation rate for different ages and different types 
of disk galaxies. The results are shown for ages of 0.1 (solid 
line), 1 (solid thickline), 3 ( dashed line) , 6 (dashed thickline) 
and 10 Gyr (dotted line). 
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2 4 6 8 
WGyrs) 



2 4 6 8 

tiow(Gyrs) 



Figure 2. The disk birthrate parameter b, the ratio of current SFR to the average past SFR, plotted against the free parameter ti ow , 
for different types of disk galaxies. The dotted line indicates the median values of birthrate parameter while the dashed lines show the 
quartiles of birthrate parameter presented by Kennicutt et al. For clarification of the plot, we define a new type of disk galaxy Sc/d. The 
data for Sc/d type come from the average between the Sc model data set and the Sd model data set. 
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Figure 3. The evolution of the gas surface densities for different ages and different types of disk galaxies. The results are shown for ages 
(top - bottom) of 0, 0.1, 1, 3, 6 and 10 Gyr (dashed line). 



2.12.5 Metallicity Profiles 

Figure [7] shows the evolution of the radial distribution of 
gas metallicity for different ages and different type of disk 
galaxies. The results are shown for ages of 0.1 (solid line) 
and 10 Gyr (dotted line). 

Note the very high value of metallicity for Sa type galax- 



ies today. The model predicts that metal-rich winds are 
needed especially for Sa type galaxies in order to produce 
reasonable values for metallicity today. This prediction is in 
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Figure 4. The evolution of stellar surface densities for different ages and different types of disk galaxies. The results are shown for ages 
(top - bottom) of 0, 0.1, 1, 3, 6 and 10 Gyr (dashed line). 




Figure 5. The radial distribution of the gas velocity at 10 Gyr for different types of disk galaxies, Sa (solid line), Sb (thick solid line), 
Sc (dashed line), Sd (thick dashed line) and Sm (dotted line). 



agreement with suggestions presented in Dalcanton (2006) 
(for more detailed comments see Section [qT2^R . 



3 Note that metallicity today does decline slightly close to the 
center. This is due to the radial profile we choose for the infall 
rate since we add most of the fresh low metallicity gas to the 
central regions. 



2.12.6 Infall Rate Profiles 



The evolution of the radial distribution of the infall rate (in 
units of M0pc -2 Gyr -1 ) is shown in Figure [8] The results 
are shown for different ages. As already mentioned the time 
when the infall switches on is defined by the free parameter 
tiow For the results presented in figure 5, we adopt for ti ow 
a value that reproduces the median value of birth parameter 
b for each type of disk galaxy. For Sa galaxies the results are 
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Figure 7. The evolution of the radial distribution of gas metallicity for different ages and different types of disk galaxies. The results 
are shown for ages of 0.1 (solid line)and 10 Gyr (dotted line). 



shown for 3 (solid line) and 6 Gyr (dashed line). Infall rates 
for Sb, Sc, Sd, and Sm galaxies are presented for ages of 
8.5 (solid line) and 10 Gyr (dashed line). Note the different 
range in the y axis for different types of disk galaxies. 

2.12.7 Disk Global Properties 

Figure [9] shows the net star formation rates (upper panel) 
and the evolution of gas fraction and total disk mass, stars 



plus gas (middle panels), for different types of disk galaxies 
respectively. Note the very high values for the infall rate for 
the Sb, Sc, Sd and Sm galaxies. Since these galaxies have 
smaller mass and the infall is switched on much later in 
comparison with the Sa type, we need high infall rates in 
order to reproduce the observed total mass density in the 
solar neighbourhood. 
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Figure 8. The evolution of the radial distribution of the infall rate for different ages and different types of disk galaxies. For Sa galaxies 
the results are shown for 3 (solid line) and 6 Gyr (dashed line). Infall rates for Sb, Sc, Sd, and Sm galaxies are presented for ages of 8.5 
(solid line) and 10 Gyr (dashed line) 




4-1 9 6-1 9 

t(yrs) 

Figure 9. Upper panel: time evolution of net star formation rate. Middle panels: time evolution of gas fraction and total disk mass 
(stars plus gas). Bottom panel: ratio of star formation rate at r = 1 kpc over the mean star formation rate. The results are shown for 
different types of disk galaxies, Sa (solid line), Sb (solid thick line), Sc (dashed line), Sd (dashed thick line) and Sm (dotted line). 



2.13 Star formation rate histories 



2.13.1 Galactocentric radius r = 1 kpc 



In this section, we present results for the star formation rate 
history at galactocentric radius r = lkpc for different types 
of disk galaxies. 



In figure(9] we show the star formation rate history at galac- 
tocentric radiusr equal to r = 1 kpc (bottom panel) for dif- 
ferent types of disk galaxies. Note, that galaxy types which 
are characterised by large infall rates (Sb, Sc, Sd, Sm) show a 
peak at the onset of infall. At time t — ti ow , when the infall 
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Table 5. Definition of numerical type T 

switches on, the accreted gas increases the local star forma- 
tion rate by a bigger factor in comparison with the increase 
in the mean star formation rate across the disk and causes 
the peak in Figure [9] This phenomenon is much stronger at 
smaller radii because of the inner 'hole' that has developed 
in the surface density of the gas at late times (see figure 
At late times, the star formation in the inner parts of the 
disk is stopped due to lack of gas, but the fresh gas comes in 
at t — ti ow to refuel the central regions and increase the lo- 
cal star formation by a big factor. For Sa galaxies, the infall 
rate is very small so the curve in figure [9] is smooth. 



3 REFINEMENTS OF THE MODEL 

The model described in the previous section contain a num- 
ber of simplifications, which we will attempt to refine in 
this section. We introduce some simple improvements to the 
outflow model (Section 13. 1[) and to the chemical evolution 
model (Section |3.2[) . 

3.1 Outflow model 

In Section \2. 81 we assumed that the k parameter is constant. 
Instead of keeping the k parameter in the outflow model con- 
stant for the different types of disk galaxies, we can assume 
that k is proportional to the ratio of bulge-to-disk circular 
velocities. This is a reasonable assumption especially if wc 
believe that AGN feedback contributes to the wind mecha- 
nism. For systems with constant density, the bulge-to-disk 
velocity ratio can be reduced to the following expression: 

^oc(^) 1/a (25) 
Vd r d 

where r e is the effective radius of the bulge and Td is the 
the disk scale length. Following MacArthur et al. (2003), 
we can use the following expression for the r J e - ratio : 

— = 0.20 - 0.013(T - 5) , (26) 

where T defines the galaxy type. The relationship between 
T and galaxy type is shown in Table [5] Note that eq. [26] 
is valid only for galaxies with numerical type T between 1 
and 7. Initial investigation of equation 1261 using a = 1 and 
a — 2, shows that our results are insensitive to this refine- 
ment of the model. We plan to investigate the implications 
of equation (|26|) in more detail in a future paper. 

3.2 Chemical evolution 

In Section \2. 101 we assumed that gas ejected from the disk 
has the same metallicity as the ISM at the time that the gas 
was ejected. Due to this assumption, the predicted metal- 
licty curve for Sa and Sb galaxies today is very high (see Fig- 
ure [7] top panel). This problem does not appear in smaller 



galaxies (Sc, Sd, Sm) because the infall rate is higher than 
the rate of star formation, thus allowing the accreted metal- 
poor gas to dilute the ISM faster than it can be enriched by 
evolving stars. A comparison with observational data (Pi- 
lyugin et al. (2004) ; Kewley et al. (2005)) leaves metal-rich 
outflows as the only viable mechanism for producing the low 
effective yields observed in gas-rich galaxies. Following Dal- 
canton (2006), we parameterise the metallicity Zsn of the 
SN ejecta as a multiple rj of the nucleosynthetic yield. For a 
Salpeter IMF -q — 6.2 — 7.1. For this paper, we adopt a value 
rj — 4.25 for Sa type galaxies and a smaller value rj — 2 
for the other galaxy types. Furthermore, in this improved 
model for chemical evolution for the accreted gas, we adopt 
a metallicity of Zf = 0.1 xZq. The equation of galactic 
chemical evolution is : 

Y, g dZ = pipdt + (Z F - Z)f - Z(r) - 1)M W (27) 

In figure \W\ we show the evolution of the radial dis- 
tribution of gas metallicity for different ages and different 
type of disk galaxies. The results are shown for ages of 0.1 
(solid line) and 10 Gyr (dotted line). Comparison between 
the model predictions and the observational data is shown in 
figure [T2l The results come from Kewley et al. (2005) (Ta- 
ble 2 in this paper). For the oxygen abundance in the solar 
neighbourhood, we adopt a value 12 + 0/H = 8.5 (Pilyugin 
et al. (2002)). The model predictions for global metallic- 
ities are calculated assuming that the disk is truncated at 
r/r d ta 5. 

Figure [TT] shows the radial metallicity distribution of 
the interstellar medium for each of our disk models after 10 
Gyr (dashed line). For the Sa model which is similar to the 
Milky Way the metallicity gradient is about dlog(Z)/dr= - 
0.055 dexkpc -1 at the solar radius. Observational values for 
the metallicity gradient of light elements X in the Milky 
Way are in the range of -0.04 <dlog(X/H)/dr< -0.08 
dexkpc -1 (Chiappini et al. (2001)). Our models predict 
that the metallicity gradient at the solar radius of the disks 
depends on the type, and increases as we go from early spi- 
rals to later types in agreement with observations (Marquez 
et al. (2002)). Our predictions for the metallicity gradient of 
later type disk galaxies are in reasonable agreement with ob- 
servations (Vila-Costas & Edmunds (1992)) although there 
is a considerable scatter in the published values. 



4 DISCUSSION 

We have presented a model of global star formation incorpo- 
rating supernova feedback, gas accretion and enriched out- 
flows in disks modelled by a multiphase interstellar medium 
in a fixed gravitational potential. 

A key prediction of this model is that star formation 
histories of different types of disk galaxies can be explained 
in a simple sequence of models which are primarily regu- 
lated by the cold gas accretion history. The distributions of 
disk birth parameters presented in Kennicutt et al. (1994) 
are reproduced using the parameter t\ ow which varies with 
the type of disk galaxy. Sa galaxies are characterised by 
quiescent evolution and a small value for the ti ow parame- 
ter whereas Sb, Sc, Sd and Sm galaxies are characterised by 
starbursts and relatively large values for the ti ow parameter. 

A description of disk evolution in which the onset and 
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Figure 10. The evolution of the radial distribution of metallicity for different ages and different types of disk galaxies. The results are 
shown for ages of 0.1 (solid line) and 10 Gyr (dashed line). 
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Figure 11. Radial metallicity distribution of the gas after 10 Gyr (dashed line). The solid lines indicate metallicity slopes at the solar 
radius. Sa (dlog(Z)/dr= -0.055 dex kpc" 1 ), Sb (dlog(Z)/dr= -0.06 dcx kpc" 1 ), Sc (dlog(Z)/dr= -0.065 dex kpc" 1 ), Sd (dlog(Z)/dr= 
-0.1 dcx kpc- 1 ) and Sm (dlog(Z)/dr= -0.11 dcx kpc" 1 ). 



duration of the gas infall history are found to be the con- 
trolling parameters is in qualitative agreement with stan- 
dard AC DM cosmology which predicts that protodisks re- 
side in dark halos with masses ~ 10 12 Mq and are in a phase 
of strong gas accretion with values >100 M0/yr (Burk- 
ert & Lin (2006)). 

The mass assembly of galaxies occurs through two main 
processes: hierarchical merging of smaller entities, and more 



diffuse gas accretion. The relative importance of the two pro- 
cesses cannot be easily found by cosmological simulations, 
since many physical parameters such as gas dissipation, star 
formation and feedback, are still unknown. 

There are at least some twenty examples of galax- 
ies which in HI show either signs of interactions and/or 
have small companions (Sancisi (1999)). This suggests that 
galaxies often are in an environment where material for ac- 
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Figure 12. The integrated metallicity is plotted versus morpho- 
logical type. Observational dat&(diamond points) come from Kew- 
ley et al. 2005, plus points correspond to the improved model 
predictions. 



cretion is available. Characteristic examples are the com- 
panions NGC 4565-4565A (Rupen (1991)) and NGC 4027- 
4027A (Phookun et al. (1992)). These companions have 
systematic velocities close to those of the main galaxy and 
HI masses less than 10% of the main galaxy. The HI picture 
suggests the capture of a gas rich dwarf by a massive system 
probably to be followed by tidal disruption and accretion of 
the dwarf. Such examples have been seen also in the Milky 
Way. The discovery of the Sagittarius dwarf galaxy (Ibata 
et al. (1994)) shows that accretion is still taking place at 
the present time. 

The model presented here predicts that metal-rich 
winds are needed especially for Sa type galaxies in order 
to produce reasonable values for the metallicity today. This 
prediction is in agreement with suggestions presented in Dal- 
canton (2006) (for more detailed comments see Section[32]). 

The major shortcoming of the present model is the fail- 
ure of the chemical enrichment model to reproduce the ob- 
served values for nuclear metallicities. Our model predicts 
very high values for metallicities close to the disk centre for 
Sa and Sb galaxies whereas the predicted nuclear metallici- 
ties for Sd and Sm galaxies are quite low. This is probably 
due to the high star formation rates that the model predicts 
in the first 10 8 yr, due to our implicit assumption that the 
entire gas disk has formed instantaneously at t = 0. The 
extremely high initial star formation rate result in highly en- 
riched central regions, so enriched that the fresh low metal- 
licity accreted gas cannot dilute the gas efficiently. This is 
the case especially in galaxies like Sa and Sb characterised 
by small infall rates. In galaxy types like Sd and Sm, the 
very high infall rates in the central regions (due to the func- 
tional form we assume for the infall model: see Sec 12.61) . 
are enough to dilute metals and result in low values for nu- 
clear metallicities. However the chemical evolution model is 
in good agreement with observations in terms of the inte- 
grated metallicities. 

In a future paper, we will develop this model further. 
Along with the improvements mentioned above, we intend 
to examine the effect of AGN feedback and compare with 
numerical simulations. 



Helly J. O, Frenk C. S., 
2006, MNRAS, 370, 645 
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